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Abstract. We implement the effective field theory for gravitating spinning objects in the 
post-Newtonian scheme at the next-to-next-to-leading order level to derive the gravitational 
spin-orbit interaction potential at the third and a half post-Newtonian order for rapidly ro¬ 
tating compact objects. Erom the next-to-next-to-leading order interaction potential, which 
we obtain here in a Lagrangian form for the first time, we derive straightforwardly the corre¬ 
sponding Hamiltonian. The spin-orbit sector constitutes the most elaborate spin dependent 
sector at each order, and accordingly we encounter a proliferation of the relevant Feynman 
diagrams, and a significant increase of the computational complexity. We present in detail 
the evaluation of the interaction potential, going over all contributing Feynman diagrams. 
The computation is carried out in terms of the “nonrelativistic gravitational” fields, which 
are advantageous also in spin dependent sectors, together with the various gauge choices 
included in the effective field theory for gravitating spinning objects, which also optimize the 
calculation. In addition, we automatize the effective field theory computations, and carry 
out the automated computations in parallel. Such automated effective field theory computa¬ 
tions would be most useful to obtain higher order post-Newtonian corrections. We compare 
our Hamiltonian to the ADM Hamiltonian, and arrive at a complete agreement between 
the ADM and effective field theory results. Finally, we provide Hamiltonians in the center 
of mass frame, and complete gauge invariant relations among the binding energy, angular 
momentum, and orbital frequency of an inspiralling binary with generic compact spinning 
components to third and a half post-Newtonian order. The derivation presented here is es¬ 
sential to obtain further higher order post-Newtonian corrections, and to reach the accuracy 
level required for the successful detection of gravitational radiation. 
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1 Introduction 

In light of the upcoming operation of second-generation ground-based interferometers world¬ 
wide, such as Advanced LIGO in the US [1], Advanced Virgo in Europe [2], and KAGRA 
in Japan [3], we may witness a direct detection of gravitational waves (GW) within the end 
of the decade, which will open a new era of observational gravitational wave astronomy. 
Inspiralling binaries of compact objects are of the most promising sources in the accessible 
frequency band of these experiments, where they can be treated analytically within the post- 
Newtonian (PN) approximation of General Relativity [4]. As the search for the GW signals 
employs the matched filtering technique, there is a pressing need to obtain accurate waveform 
templates, using the Effective-One-Body (EOB) formulation to model the continuous signal 
[5]. 

Recently the fourth PN (4PN) order correction has been completed for the non spinning 
case [6], and it is necessary to reach a similar accuracy level for the spin dependent case, 
as such compact objects are expected to be rapidly rotating. Of the spin dependent sectors 
the spin-orbit sector contributes the leading order (LO) spin-dependent PN correction, and 
represents the most dominant spin effects. The LO spin-orbit correction at the 1.5PN order 
has been obtained by Tulczyjew already at 1959 [7], see erratum therein, and later in the form 
of a higher-order Lagrangian used here by Damour [8]. Yet, the next-to-leading order (NLO) 
spin-orbit interaction has been approached much later, first at the level of the equations of 
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motion (EOM) in [9], and [10], and then within the ADM Hamiltonian formalism [11]. The 
next-to-next-to-leading order (NNLO) spin-orbit interaction was first derived in Hamiltonian 
form in [12, 13], building on [14, 15], and then at the level of the EOM in [16, 17]. 

In this work we apply the effective field theory (EFT) for gravitating spinning objects 
in the PN scheme [18] at the NNLO level, which was already tackled in the spin dependent 
sector using EFT techniques in [19]. We derive the NNLO gravitational spin-orbit interaction 
potential at the 3.5PN order for rapidly rotating compact objects. The EFT for gravitating 
spinning objects [18] builds on the novel, self-contained EFT approach for the binary inspiral, 
which was introduced in [20, 21]. The seminal works in [22, 23] already treated spin in an 
action approach in flat and curved spacetime, respectively. Then an extension to spinning 
objects of EFT techniques was approached in [24], which later adopted a Routhian approach 
from [25], and leaves the temporal components of the spin tensor in the final results. Yet, 
these components depend on field modes at the orbital scale, and they must be eliminated in 
order to obtain physical observables. The EFT for gravitating spinning objects [18] actually 
obtains an effective action at the orbital scale by integrating out all the field modes below this 
scale. Moreover, it actually enables the relation to physical observables: the physical EOM 
are directly obtained via a proper variation of the effective action [18, 26]. Furthermore, it 
also enables to obtain the corresponding Hamiltonians in a straightforward manner from the 
potentials derived via this formulation [18]. Indeed, from the potential, which we obtain here 
in a Lagrangian form, we derive the corresponding NNLO spin-orbit Hamiltonian, and then 
compare our result to the ADM Hamiltonian in [12, 13]. We arrive at a complete agreement 
between the ADM and EFT results. 

The spin-orbit sector constitutes the most elaborate spin dependent sector at each order, 
see [18] for the LO and NLO levels, and [19, 27] for the other sectors at the NNLO level. 
Accordingly, we encounter here a proliferation of the relevant Feynman diagrams, where there 
are 132 diagrams contributing to this sector, and a significant increase of the computational 
complexity, e.g. there are 32 two-loop diagrams here. We also recall that as the spin is 
derivative-coupled, higher-order tensor expressions are required for all integrals involved in 
the calculations, compared to the non spinning case. Yet, the computation is carried out 
in terms of the “nonrelativistic gravitational” (NRG) fields [28, 29], which are advantageous 
also in spin dependent sectors, as was first shown in [30], and later also in [18, 19, 31, 32]. We 
also apply the various gauge choices included in the EFT for gravitating spinning objects [18], 
which also optimize the calculation. In addition, we automatize the EFT computations here, 
and carry out the automated computations in parallel, where we have used the suite of free 
packages xAct with the Mathematica software [33, 34]. Such automated EFT computations 
would be most useful to obtain higher order PN corrections. It should be stressed that in 
order to obtain further higher order results, all lower order results are required, consistently 
within one formalism, and so also in that respect the derivation presented in this work is 
essential. Finally, we provide Hamiltonians in the center of mass frame, and complete gauge 
invariant relations among the binding energy, angular momentum, and orbital frequency of 
an inspiralling binary with generic compact spinning components to 3.5PN order. 

The outline of the paper is as follows. In section 2 we briefly review the EFT for 
gravitating spinning objects in the PN scheme, and present the Feynman rules required 
for the EFT computation. In section 3 we present the evaluation of the NNLO spin-orbit 
interaction potential, going over all contributing Feynman diagrams, and giving the value 
of each diagram. In section 4 we present the NNLO spin-orbit potential EFT result, and 
from it we obtain the corresponding EFT Hamiltonian. We then compare our result to the 
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ADM Hamiltonian, where we resolve the difference between the Hamiltonians, using higher 
order PN canonical transformations, and arrive at a complete agreement between the ADM 
and EFT results. We also present all relevant Hamiltonians to 3.5PN order in the center 
of mass frame. In section 5 we provide the complete gauge invariant relations among the 
binding energy, angular momentum, and orbital frequency of an inspiralling binary with 
generic compact spinning components to 3.5PN order. In section 6 we summarize our main 
conclusions. Finally, in appendix A we provide the additional irreducible two-loop tensor 
integrals required for this work. 

2 The EFT for gravitating spinning objects in the PN scheme 

In this section we present the effective action, and the Feynman rules, which are derived 
from it, and are required for the EFT computation of the NNLO spin-orbit interaction. We 
employ here the “NRG” fields, as applied with spin in [18, 19, 30-32]. Here, we briefly review 
and build on [18, 19, 26, 30, 31], following similar notations and conventions as those that 
were used there. Hence we use c = 1, r]^i, = Diag[l, —1, —1, —1], and the convention for the 
Riemann tensor is ~ scalar triple product 

appears here with no brackets, i.e. a x b ■ c = {a x b) ■ c. The notation = f is used 
for abbreviation. In fact, the generic d-dimensional dependence can be and is suppressed in 
what follows, and d can be set to 3, except for computations which involve loops, where only 
the d dependence from the generic d-dimensional Feynman integrals, see appendix A in [19], 
should be considered. 

First, in terms of the “NRG” fields the metric reads 
_ ( \ 

_ / l + 2^ + 2^^-l- 1^3 -Aj - 2Aj(j) - 2Aj(j)‘^ \ 

y —Ai — 2Ai(f) — 2Aj(l>^ — 6ij -\- 2<p5ij — cTij — 2(jP‘5ij -\- 2(f>(jij -|- AiAj -|- J 

( 2 . 1 ) 

where we have written the approximate metric in the weak-field limit up to the orders in the 
fields that are required for this sector. 

We recall that the effective action, describing the binary system, is given by 

2 

S = S8 + E^(I)pp" P'2) 

1 = 1 

where Sg is the pure gravitational action, and is the worldline point particle action for 

each of the two particles in the binary. The gravitational action is the usual Einstein-Hilbert 
action plus a gauge-fixing term, which we choose as the fully harmonic gauge, such that we 
have 

= 5eh + 5gf = I d^x^R+ ^ I (2.3) 

where T^ = 

From the gravitational action we derive the propagators, and the self-gravitational ver¬ 
tices. The “NRG” scalar, vector, and tensor field propagators in the harmonic gauge are 
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then given by 


r Ak-{x\-X2) 

- ={4>{xi) (f){x 2 ))= 47rG J -p- S{ti-t 2 ), (2.4) 

r ik-(xi-X2) 

- =(^Ai{xi) Aj{x2)) = -IGttG 5ij -p-(5(ti -t2), (2.5) 

p _^ 2 ) 

={(7ij{xi)akiix2)) = 32TrG Pij-ki —p- S{ti-t 2 ), (2.6) 

where Pij-^ki — 2 P ^u^jk ‘^SijSki)• 

The Feynman rules for the propagator correction vertices are given by 

-^- =8^ / 

---X —- = -^ I d^x idtA,)\ (2.8) 

X = 12 ^ / - {dt(Tiif] , (2.9) 


where the crosses represent the self-gravitational quadratic vertices, which contain two time 
derivatives. 

The Feynman rules for the three-graviton vertices required for the NNLO of the spin- 
orbit interaction are given by 






1 


J d'^x (j) {diAj {diAj 


djA,) + (d^Ai)^) , 


( 2 . 10 ) 


1 

GdvrG 



2(7jj [diAj^djA]^ -\- d]^Aid]fAj 


‘^‘dkAidjAjf + 2diAjdkAf^) 


-Gkk (diAj {diAj - djAi) + {diAif'^ , 


( 2 . 11 ) 


IGttG J 

^ d^x 

{2aijdicl)dj(i) - ajjdicjrdict)) + (an {dtckf - lGdtaii(t>dt(t^ 

(2.12) 

1 

47rG 

/a 

{Aidi(j)dt(j)) , 

(2.13) 

00 

II 

d'^x [2aij {di(j)dtAj - dt4>diAj) - ajj {di(f)dtAi - dt(t)diAi)] , 

(2.14) 

= ' / 
IGttG J 

d^x {AidiAjdtAj) , 

(2.15) 

1 

27rG 


(t>{dt4>f , 

(2.16) 


where the first two vertices are stationary, and can be read off from the stationary Kaluza- 
Klein part of the gravitational action. The next five vertices are time dependent, and contain 
up to two time derivatives. 
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The Feynman rule for the four-graviton vertex required to the order considered here is 
given by 


where this vertex is stationary. 

Next, we recall that the minimal coupling part of the point particle action of each of 
the particles with spins is given by 



5pp = / dX 


/' 


—mV w 


2 ^ 


Dp^ 
p2 DX 


(2.18) 


which is covariant, as well as invariant under rotational variables gauge transformations [18]. 
Here A is the affine parameter, = dx^/dX is the 4-velocity, and are the angular 

velocity and spin tensors of the particle, respectively [18]. We parametrize the worldline 
using the coordinate time t = i.e. X = t, so that we have for = dx^/dX: = 1, 

u® = dx'^/dt = u®. Since the spin-orbit interaction is linear in the spins, only the minimal 
coupling part of the action, i.e. that which appears in eq. (2.18) is required. We stress that 
in the spin-orbit sector both mass and spin couplings play central roles in the interaction. 

Let us then present first the mass couplings required for this sector. The Feynman 
rules of the one-graviton couplings to the worldline mass required for the NNLO spin-orbit 
interaction are given by 


/' 


= — m / dt 4> 


- = m / dt AjV 


/' 
4»/ 


. 3 2 74 

l + r +' 


. 1 2 34 

1 + + 


dt aijV^v^ 


'^ + ^v + 


(2.19) 

( 2 . 20 ) 

( 2 . 21 ) 


where the heavy solid lines represent the worldlines, and the spherical black blobs represent 
the masses on the worldline. The ellipsis denotes higher orders in u, beyond the order 
considered here. 

For the two-graviton couplings to the worldline mass required here, we have the following 
Feynman rules: 
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(2.24) 



dt (paijV^'v^ 


[1 + ...]. 


Finally, for the three-graviton couplings to the worldline mass required here, we have 
the following Feynman rules: 





(2.25) 



dt (j)^Aiv'^ [1 -f • • • ]. 


(2.26) 


Let us go on to the spin couplings required for this sector. These are given here in 
terms of the physical spatial components of the local spin variable in the canonical gauge 
[18]. First, we have contributions from kinematic terms involving spin without field coupling 
[18]. To the order we are considering, these are given by 

Lkin = —S ■ Q + —S • F X a ^1 + —> (2.27) 

where Sij = eijkSk, = eijk^k, ^ijk is the 3-dimensional Levi-Civita symbol, and a* = uL 
We recall that all indices are Euclidean. 

The required Feynman rules of the one-graviton couplings to the worldline spin are thus 


\-= / dt 


(^ijkSk ( i^diAj + {diAj - djAi) ( 3 j v^dtAj ( 1 + 


-|- Aiv'' 


j dt eijkSkV^ ^dj(f) + - a^cj) (2 -k 3u^)^ 


(2.28) 


(2.29) 


/ 


= dt 


^^ijkSk (diaji v'‘ +i ]^diaira - diaimj v^v^v'^ - ^dtau v^v'- 


- an I v^a — V 


(2.30) 


where the (gray) oval blobs represent the spins on the worldlines. 

For the two-graviton couplings to the worldline spin, the Feynman rules required here 

are: 



= dt [2eijkSk {diAj (j) -|- Aj dt(j)v'‘)] 


(2.31) 
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dt 


^(^ijkSk^il {djA-i dlAj') 


dt 


2^ijkSk (^djAi Ai + Ai dtAjJ 


(2.33) 



= J dt [4:eijkSk<p [di(f)v^ + v'^a\ 



= dt 


(-ijkSkfyu (djfpv’- + dicpv^ 


(2.34) 


(2.35) 


Finally, we also have to include three-graviton spin couplings. For three-graviton cou¬ 
plings to the worldline spin, the only Feynman rule required here is 



= j dt [AeijkSkdiAj (j)^ 


(2.36) 


Note that similarly to the NLO, where the two-scalar spin coupling is absent, here at NNLO 
the three-scalar spin coupling vanishes due to the use of the “NRG” helds, and our gauge of 
the rotational spin variables. Its appearance is then deferred to higher PN orders, which is 
advantageous. 


3 Next-to-next-to-leading order spin-orbit interaction 

In this section we evaluate the relevant two-body effective action by its diagrammatic expan¬ 
sion. As explained in [19], in the NNLO spin-orbit potential, which is evaluated at 3.5PN 
order, we have diagram contributions up to order coming from all 12 possible topologies 
appearing at these orders, as displayed in figures 1-6 below: one topology at 0{G), two at 
O(G^), and nine topologies at 0{G^). For the construction of the Feynman diagrams we use 
the Feynman rules presented in the previous section, see [19] for more detail. We have a total 
of 132 diagrams contributing to the NNLO spin-orbit interaction. Eventually, we will also 
have contribution from terms with higher order time derivatives coming from the point-mass 
2PN order [35], and up to NLO spin-orbit [18] sectors. 

We denote r = xi — X2, r = |f), and ft = ^. Labels 1 and 2 are used for the left and 
right worldlines in the figures, respectively. All of the diagrams should be included together 
with their mirror images. Accordingly, the (1 •H- 2) notation stands for a term, whose value 
is obtained under the interchange of particles labels. Finally, a multiplicative factor of f dt 
is omitted from all diagram values. 

3.1 One-graviton exchange 

For the NNLO spin-orbit interaction we have 8 one-graviton exchange diagrams as shown in 
figure 1. In addition to the one-graviton exchange diagrams, which already appeared in the 
NLO spin-orbit sector [18, 31], new diagrams are added here by inserting further propagator 
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Figure 1. NNLO spin-orbit Feynman diagrams of one-graviton exchange. 


correction vertices. Tensor Fourier integrals of up to order 5 are required here, due to the 
derivative-coupling of spin, which makes the computations heavier, see [19], and appendix A 
there. 

At the NNLO level we inevitably obtain terms with higher order time derivative, i.e. with 
accelerations and precessions, and these are all kept until they are treated rigorously in the 
resulting action [18, 26], see section 4 below. Finally, we recall that there are several ways 
to evaluate diagrams with time derivatives, which differ only by total time derivatives. Our 
convention for their evaluation is, that time derivatives from the spin couplings are taken on 
their respective worldlines, whereas those from the propagator correction vertices are taken 
symmetrically on the worldlines. 

The values of the one-graviton exchange diagrams are given as follows: 
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Fig. 1(c) 


Fig. 1(d) 


Fig. 1(e) 
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Figure 2. NNLO spin-orbit Feynman diagrams of two-graviton exchange. 


— 5i X n • a 2 Vi ■ V 2 ) (3-8) 

3.2 Two-graviton exchange and cubic self-interaction 
3.2.1 Two-graviton exchange 

For the NNLO spin-orbit interaction we have 15 two-graviton exchange diagrams, as shown 
in figure 2 , where they either contain a two-graviton spin or mass coupling. 

The two-graviton exchange diagrams also require only tensor Fourier integrals. We 
encounter here two-graviton exchange diagrams, which involve time derivatives, either from 
the spin couplings or from propagator correction vertices. 

The values of the two-graviton exchange diagrams are given in the following: 
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3.2.2 Cubic self-interaction 


For the NNLO spin-orbit interaction we have 49 cubic self-interaction diagrams, as shown in 
figure 3 , where the cubic vertices contain up to two time derivatives. 

The cubic self-gravitational interaction diagrams require first the application of one-loop 
tensor integrals up to order 3 , in addition to the Fourier tensor integrals, see appendix A in 
[ 19 ]. Here, we encounter time derivatives from the spin couplings, the propagator correction 
vertices, and the time dependent cubic self-gravitational vertices. 

The values of the cubic self-interaction diagrams are given as follows: 
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Figure 3. NNLO spin-orbit Feynman diagrams of cubic self-interaction. 
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2G‘^ml ^ 


5i X n • F 2 ( 4 Fi • V 2 + 3?;| — I 6 F 1 • fiv 2 ■ n — 16(F2 • n)^) 


+ 45 i X Fi • F2F2 • fi 


1 SG'^rrin 

r ^ ^ 1 

+ 

^ J.2 

25i X fi ■ 02 F 2 • n — 5i X F 2 • 02 


8G2 


m? 


V 2 ■ fiSi X fi ■ V 2 


(3.69) 


2G^mim2 


5i X n • Fi ( 3 F 1 • V2 — 4Fi • fiv2 ■ n) — 35i x n • F 2 {v\ + Fi • F 2 
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Figure 4. NNLO spin-orbit 


Feynman diagrams of order with no loops. 


— Avi ■ nv 2 ■ n — 4{vi ■ n)^) — 25i x vi ■ V 2 [Svi ■ n + V 2 ■ n) 
2G^mim2 


+ 


Fig. 3(gl) = - 


- Si X vi ■ a2 

+ 35i X vi ■ V2 
2G‘^m? ■■ 


3S'i X n ■ aiV 2 ■ n + 3S'i x n ■ 02 F 1 • n -|- S*! x oi • 172 
2G‘^mim2 


+ 


5i X n • 'i;ii ;2 • n — 45i x n • V 2 V 1 ■ n 


(3.70) 


S'! X n . iTi {4vi ■ V 2 — V 2 — 16i;i • nv2 ■ n -|- 2{v2 ■ n)^) 


— Si X vi ■ V 2 (4Fi ■ n — V 2 ■ n) 


■■ 


4G2 


m 


2 r 


2S'i X n ■ V 1 V 2 ■ n + Si x vi ■ V 2 


2Si X fi ■ aiV2 ■ fi + Si X ai ■ V2 

(3.71) 


T7- o^ 2G2mim2 
Fig. 3(g2) =- -3 - 


-b 4S'i X vi ■ V 2 V 1 ■ fi 


Si X n ■ vi [vf — 4i;i • V 2 + IQvi ■ nv2 ■ n — 2{vi ■ n)^) 
2G^mim2 


+ 


— 25i X Si ■ V2 + 


3.3 Cubic in G interaction 


2G^mim2 


Si X fi ■ ai{yi ■ fi — 4 v 2 ■ n) — Si x vi ■ di 

Si X fi ■ vi{yi ■ fi — 4 v 2 ■ n) — 2Si x vi ■ V 2 

(3.72) 


3.3.1 Three-graviton exchange 

For the NNLO spin-orbit interaction we have 7 diagrams at order G^ with no loops, as shown 
in figure 4. These three-graviton exchange diagrams are constructed with either one-, two-. 


or three-graviton spin couplings. 

The values of these diagrams are given by 

Fig. 4(al) = - 8 *^ Si ■ V 2 x n, (3.73) 

Fig. 4(a2) = - 8 ^ Si ■ V 2 x n, (3.74) 

Fig. 4(a3) = - 2 ^ Si • F 2 x n, (3.75) 

Fig. 4(a4) =2 ^ Si • Fi x n, (3.76) 

Fig. 4(bl) = - • F 2 X n, (3.77) 
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Figure 5. NNLO spin-orbit Feynman diagrams of order with one loop. 


Fig. 4(b2) = - ^ ■V2xn, (3.78) 

Fig. 4(b3) X n. (3.79) 

3.3.2 Cubic self-interaction with two-graviton exchange 

For the NNLO spin-orbit interaction we have 21 diagrams at order with one loop, as shown 
in figure 5. These diagrams contain both cubic self-interaction and two-graviton worldline 
couplings. 

The values of these diagrams are given by 


- 

Fig. 5(al) =32 • i ;2 x n, 

(3.80) 

F.g. 5(a2) 

(3.81) 

Fig. 5(a3) 

(3.82) 

F.g. 5(a4) = 8 ®’‘”r"Sl-- 2 Xn. 

(3.83) 

Fig. 5(a5) =2^ Si ■ V 2 x n, 

(3.84) 

F.g. 5(a6) 

(3.85) 

Fig.5(a7)- 8 ®'‘T’""-Si'”iX". 

(3.86) 
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Fig.5(a8)- Si ■*><«. 

(3.87) 

Fig.5(a9)- 

(3.88) 

(~T^Tn^ —k 

Fig. 5(bl) = • U 2 X n. 

(3.89) 

Fig. 5(b2) - ‘^Si-vixn, 

(3.90) 

Fig. 5(b3) - ^ 

(3.91) 

Fig. 5(b4) =0, 

(3.92) 

Fig, 5{b5) 

(3.93) 

Fig.5(b6)=T“’"‘’"L-,.ftxfl, 

2 

(3.94) 

Fig,5{b7)=T“"‘|’"^S,.F,gjI, 

2 r* 

(3.95) 

Fig. 5(cl) - 2 • U 2 X n. 

(3.96) 

F.g. 

(3.97) 

Fig.5(c3)= 2°“’"r^S..ax3. 

(3.98) 

Fig. 5(c4) - ^^i-uixn, 

2 

(3.99) 

Fig. 5(c5) - 1 Si-vixn. 

2 

(3.100) 


Note that the total value of the diagram in hgure 5(b4) equals 0, although it does not stand 
for a short distance contribution. 

3.3.3 Two-loop interaction 

For the NNLO spin-orbit interaction we have 32 two-loop diagrams at order G^, as shown in 
figure 6. These diagrams contain two cubic vertices or one quartic vertex, and even include 
cubic vertices with time dependence. As explained in [19], they contain two-loop Feynman 
integrals of three kinds: Factorizable, nested, and irreducible. The factorizable two-loop 
diagrams do not contribute at the NNLO level, and they yield here purely short distance 
contributions, of the form which are contact interaction terms. For other two-loop 

diagrams calculations should be made, keeping the dimension d general, and the limit d —>• 3 
is only taken in the end. 

For the irreducible two-loop diagrams, which are the most complicated, irreducible two- 
loop tensor integrals of order 3 are encountered here. These are reduced using the integration 
by parts method to a sum of factorizable and nested two-loop integrals, as explained in [19], 
and see appendix A there. In addition to the irreducible two-loop tensor integrals, which 
were given in appendix A of [19], eqs. (All), (A12) there, two further irreducible tensor 
integrals are required here, and we provide them in appendix A below. 

The values of the two-loop diagrams are given in the following: 
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(al) (a2) 







(b7) 



(bll) (bl2) 





(bl6) 



Figure 6. NNLO spin-orbit Feynman diagrams of order with two loops. 
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Fig. 6(al) =0, 

(3.101) 

Fig. 6(a2) =0, 

(3.102) 

Fig. 6(a3) =0, 

(3.103) 

Fig. 6(a4) =0, 

(3.104) 

Fig. 6(a5) =0, 

(3.105) 

Fig. 6(bl) = - 32^-^^5i ■ V2 x ft, 

(3.106) 

Fig. 6(b2) = 5i • 112 X n, 

5 

(3.107) 

Fig. 6(b3) - 1 Si-vi xn, 

5 

(3.108) 

Fig.6(b4)=|0'’"|’"^S..axfl, 

5 

(3.109) 

Fig. 6(b5) =2‘^^7^5 i -Fi x n, 

(3.110) 

Fig. 6(b6) = Si’Vixn, 

5 

(3.111) 

Fig,6(b7)=f“i”=S:.axn, 

5 

(3.112) 

Fig. 6(b8) - ^^5i-i;2Xn, 

5 

(3.113) 

Fig. 6(b9) - 1 Si-V2Xn, 

5 

(3.114) 

Fig.6(bl0)- 

(3.115) 

Fig. 6(bll) =2^ ■V 2 X n, 

(3.116) 

Fig. 6(bl2)=-iO’i”^S..ftxS. 

5 

(3.117) 

Fig.6(bl3)=l°'"'f’"^S..axfl, 

5 

(3.118) 

9 _* 

Fig. 6(bl4) = 2 5i-i;2Xn, 

5 

(3.119) 

Fig. 6(bl5) - 2 5i • Hi X n, 

5 

(3.120) 

Fig. 6(bl6)-®'‘”]”^Si.FixS, 

(3.121) 

Fig.6(cl)= 

(3.122) 

Fig. 6(c2) =4°“™r^S,.FiXn, 

(3.123) 

Fig. 6(c3) =12°“’";’"^S..FiXn, 

(3.124) 
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Fig. 6(c4) = • i;2 X n. 

(3.125) 

Fig. 6(c5) X n 2F2 X n] , 

(3.126) 

Fig. 6(c6) Si ■ [—12i;i x n -|- 10u2 x n] , 

(3.127) 

Fig. 6(dl) =0, 

(3.128) 

Fig. 6(d2) =0, 

(3.129) 

Fig. 6(el) =16^-^Si ■ V 2 x n, 

(3.130) 

T7- Cl oi G^mlm2 ^ ^ ^ 

Fig. 6(e2) = Si • 112 X n, 

5 

(3.131) 

Fig, 

5 

(3.132) 


4 Next-to-next-to-leading order spin-orbit potential and Hamiltonian 


Summing up all of the Feynman diagrams from the previous section, we obtain the NNLO 
spin-orbit interaction potential for a binary system of compact spinning objects. We split 
the potential into several pieces according to the number and order of higher-order time 
derivatives as follows: 


.O ( 0 ) (b ( 2 ) ( 3 ) ( 4 ) 

^NNLO = 'F + 'F + F + 'F + 1/ 


(4.1) 


The ordinary part of the potential, which does not contain higher-order time derivatives, 
reads 


( 0 ) 

V=- 


Gm2 
4r2 L‘ 


—* /O OO 0 / 1/1 o 

Si X fi ■ vi [5viVi ■ V2 — ^2 + 4i;i • V 2 V 2 — 3vi — 3v2 + 6ui • nviV2 


n 


— 6ui • nv2 ■ nvi ■ V 2 + bui • nv2 ■ nv2 + 3v2{vi ■ n)^ -|- 3vl{v2 ■ n)'^ — 15(ui • n)^(u2 • n)^) 
-|- Si X n • V 2 {v‘fv 2 — 2vi ■ V 2 V 2 — 2(vi ■ U 2 )^ -|- 3v2 — 6vi ■ nv2 ■ nv 2 — 3u|(i;i • n)^ 

— 3ui(i;2 • n)^ -|- 15('Ui • n)^(u 2 • n)^) -|- Si x ui • i; 2 (i;i • nvf — 2vi ■ nvi ■ V 2 + 2vi ■ nu| 

-|- 2v 2 ■ nv2 — &vi ■ n{v2 ■ n)^) 


+ 


4^3 


Si X n • i;i(l3ui — 41i;i • V 2 + 28^1 — 16ui • nv 2 ■ n + 8{vi ■ n)^ -|- 12(i;2 • re)^) 


— 2Si X n ■ V 2 {vf -|- 7vi ■ V 2 — 8^2 — 56i;i • nv 2 ■ n — 4(i;i • n)^ -|- 62(i;2 • n)^) 

2G‘^mim2 


+ Si X vi ■ i;2(7i;i • n — 38v2 ■ n) 


Si X n • ui [2vf -j- vi ■ V2 — 3V' 


+ 12 i;i • nv 2 ■ n — 16 ('Ui • n )2 -|- 4(i;2 • n)^) — 2 Si x n ■ V 2 {vi ■ V 2 — 
-|- 4i;i • nv 2 ■ n — 4(i;i • n)^) — Si x ui • i; 2 ( 2 i;i • n — 5v2 ■ n) 


+ 


G^m\m2 


Si X n ■ vi — Si X n ■ V2 


+ 


5G^mim2 


Si X n • ui — Si X n • U 2 


+ 


G^ 


TTln 


Si X n • ui — Si X n • U 2 


+ (lf^2). 


(4.2) 
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Due to the large number of terms, it makes sense to sort the terms with a single higher-order 
time derivative into 


( 1 ) 

V=Va + Vs, 


(4.3) 


where 


^o. = - —Si X vi • div\ 

lb 


+ 


Gm2 
4r L‘ 


Si X n ■ vi (4i;i • oii;2 • n — 4v2 ■ nSi ■ V 2 + ai ■ — uf 02 • n -|- 4i;i • nvi ■ 0,2 


— Avi ■ nv 2 ■ 3,2 + 3a2 • n{vi ■ n)^ — 3ai • n{v 2 ■ u)^) + 25i x n ■ Si (vfv 2 ■ n — V 2 ■ nvi ■ V 2 
+ vi ■ nv 2 + V 2 ■ nv 2 — 3ui • n{v 2 ■ n)^) — 5i x ui • ai (l2uf — 12i;i • V 2 + 5^2 

— 4i;i • nv 2 ■ n — {v 2 ■ n)^) -|- 5i x n • 'i;2 [ 2 v 2 ■ nai ■ V 2 — 3i ■ nuf + vla 2 ■ n 

— 2 vi ■ nvi ■ 02 -|- 4 ui • nv 2 ■ 02 — 3 a 2 • n{vi ■ n)^ -|- 3 ai • n{v 2 ■ n)^) 

— 25i X ui • 'i;2 (bui • oi — oi • '02 — ui • 02 + 2v2 ■ a 2 — vi ■ no2 • n) 

— Si X Si ■ V2 (buf — 4 i;i • V2 + 5 v 2 — 6 ui • nv2 ■ n + {v2 ■ u)^) -|- 25 i x n • 02 (v‘fv2 ■ n 
+ vi ■ nv2 — 2>V2 ■ n{vi ■ n)^) -|- iSi x ui • 02 {y\ — 2v2 + 4 i;i • nv2 ■ n — {vi ■ n)^) 

— 5i X i;2 • a2{vj - {vi ■ n)^) 


2G^mim2 


Si X n ■ vi (2ai ■ n + 02 ■ n) + Si x n ■ ai (3i;i • n + 2v2 ■ u) 


— 25i X n ■ 02^1 • n + 2Si x ai ■ V 2 + 2Si x vi ■ a 2 


G'^m. 


2 r 
2 




2Si X n ■ vidi ■ n 


+ 45i X re • oi (ui • n — U 2 • fi) + 25i x n ■ V 2 (oi • re -|- 14a2 • re) — 25i x re • 02 (ui 
— 18 u2 • n) + 14S'i X ui • ai -|- 27Si x ai ■ V 2 + Si x vi ■ a 2 + 125i x V 2 ■ 32 
+ (1^2), 


re 


(4.4) 


Vs=— 

^ 2 r 


Si X n ■ vi [vfv 2 ■ n — V 2 ■ nvi ■ V 2 + vi ■ reu^ + V 2 ■ nv 2 — 3i;i • n{v 2 ■ re)^) 


— 5 i X re • i;2(ui • nv2 + V2 ■ nv2 — 3 vi ■ n{v2 ■ re)^) — S'! x ui • V2{^vl + 2^2 

— 3i;i • nv 2 ■ n + {v 2 ■ re)^) 


G‘^mim2 


bS*! X n ■ vi (3i;i ■ n — V 2 ■ n) — 2Si x n ■ V 2 (7vi ■ n — 5v2 ■ re) 


-I- 75i X ui • i;2 


G^rrej 

^l’2 


5i X re • ui(3i;i • re — 4i;2 • re) -|- 2Si x n ■ V 2 {vi ■ 


re 


-|- 15i;2 • n) + 27Si x vi ■ V 2 -|- (1 -f-)- 2). 

The piece with two higher-order time derivatives is given by 


(2) 1 
V= - -Gm2 
4 


5i X re • oi {y 2 — (u2 • n)'^) + x vi ■ 3iV2 ■ n + b5i x 3i ■ V 2 V 2 


n 


(4.5) 
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— Si X ft ■ 0,2 (^1 — (wi • n)^) — 2S'i X vi ■ 02^1 • n 
1 


--Gm2 

-\Gm2 


Si X fi ■ vi {y 2 — {v 2 ■ n)'^) — Si x n ■ V 2 {^2 ~ (^2 • n)"^) + 65i x vi ■ V 2 V 2 ■ n 
Si X fi ■ vi (5ai ■ 02 + oi ■ na 2 ■ n) + 25i x n ■ oi (2vi ■ 02 — 2v2 ■ 02 
+ vi ■ na 2 ■ n) — Si X vi ■ 0 , 10,2 ■ n — Si x n ■ V 2 (3ai ■ 02 + oi ■ 002 ■ n) 

+ 5 i X ai • 1)202 ■ n — 2Si x n ■ 02 [oi ■ V2 + oi ■ nv2 ■ n) — Si x vi ■ 0,20,1 ■ n 

— 2Si X oi ■ 0,2 (4^1 ■ fi — V2 ■ -\- Si X V2 ■ «20i • n 

— -Gm2 2 Si X n ■ vi ( 2 vi ■ 02 — 2 ,V 2 •02 + ^1- na2 ■ n) + 2 Si x n ■ oi (v2 — {v2 ■ n)'^) 

+ 4S'i X vi ■ diV 2 ■ fi — 2Si X fi ■ V 2 [vi ■ 02 — 2v2 ■ 02 + Vi ■ na 2 • n) + 2Si x vi ■ V 2 a 2 ■ n 

+ 12iS'i X oi ■ V 2 V 2 ■ n + 25i x n ■ 0,2 {y 2 — 2vi ■ fiv 2 ■ n) — 451 x 'Ui • 02 (2'Ui ■ n — V 2 ■ n) 

+ 25 i X V 2 ■ a 2 Vi ■ fi 


G‘^mim2 


Si X fi ■ oi + 65i X fi ■ 0,2 


17G2 


m 


2 . 


4r 


■5i X fi • 02 + (1 0 2). 


(4.6) 


The contribution with three higher-order time derivatives reads 


(3) 1 

F=-Gm2r 


(5i X n • 0102 • n + 5i X n • 0201 • n + 75i x ai • 02 -|- 5i x oi • 0 , 2 ) 


+ (25i X n • 02^1 • n -|- 25i x vi ■ 0 , 2 ) + (5i x n ■ vid 2 ■ n — Si x n ■ '0202 • n 
— 25i X n • d2V2 ■ n + 75i x ffi ■ 0,2 — Si x V2 ■ 02) 


-|- (25i X n • 0102 • n + 145i x oi • 02) -|- (1 2), 


(4.7) 


and finally, even a term of four higher-order time derivatives appears, and reads 

( 4 ) 1 li 

V=-Gm 2 r^Si x n ■ 02 + (1 <->■ 2). (4.8) 

It is obvious that the higher order time derivatives of the velocity and spin bloat the potential. 
These can be handled at the level of the EOM through a substitution of lower order EOM. 
However, it is often more useful to perform the elimination of higher order time derivatives 
at the level of the potential, and to also transform to a Hamiltonian. This will make the 
result considerably more compact. 

For the reduction of higher order time derivatives we follow the procedure outlined in 
[26, 36], and its explicit extension for spin variables in [26]. It should be stressed that this 
procedure is in general not equivalent to a substitution of EOM at the level of the potential, 
but rather to a redefinition of variables, which combines with the higher order time derivatives 
to a total time derivative. Yet, as long as this redefinition contributes only linearly to the 
level of approximation, its result is equivalent to an insertion of the lower order EOM, up to 
a total time derivative. This total time derivative is the same that arises for the derivation 
of equations of motion through linear variation of variables. Variation and redefinition of 
variables are essentially the same. 
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We are eliminating the higher order time derivatives successively at each PN order, 
starting with the LO spin-orbit potential at 1.5PN order. As in [7], also see [8, 37], and also 
discussed in [18], the variable shift reads 


yi^yi + 


1 

2 mi 


Si X ui, 


(4.9) 


and similar for 2/2 • This shift corresponds to an insertion of EOM into the LO spin-orbit 
potential, where the EOM was derived from the complete 3.5PN order potential. The redefi¬ 
nition in eq. (4.9) is linear in the spin, so that its square does not contribute to the spin-orbit 
sector. After this step the potential, and thus the EOM are changed. The higher order 
time derivatives at the next PN orders are then eliminated using these modified EOM (the 
corresponding redehnitions of variables are omitted here due to their length). The remaining 
higher order time derivatives start to appear at 2PN order, so that the redefinitions would 
contribute quadratically, and thus the insertion of EOM breaks down, only at the 4PN order. 
That is, this successive insertion of EOM to remove higher order time derivatives is still valid 
to linear in spin and 3.5PN order. 

It should be noted that this procedure is different from inserting the EOM in the com¬ 
plete 3.5PN order potential in a single step. This would correspond to a redefinition of the 
position, which contains among others a contribution from the LO spin-orbit and the 2PN 
order point-mass potentials. The quadratic contribution of this redefinition therefore leads 
to NNLO spin-orbit terms. That is, the method of inserting EOM breaks down in this case, 
in contrast to the successive elimination discussed above. Furthermore, both procedures can 
lead to different (yet equivalent) results, since the total time derivatives generated by the 
successive linear redefinitions and by the corresponding single quadratic redehnition are in 
general different. 

Next, we can perform a Legendre transformation to obtain a Hamiltonian. For that we 
need to replace the velocity in terms of canonical momenta, which reads 


=P\ - \p\pI + [ - ^P\ + 7P2 + n*n • # 2 ] - 

+ ^ [ - heikjP2Sin • Pi + AeikjP2Sin • p2 + QeikjU^Sln ■ pin ■ p2 

-|- 2eikjn^S 2 P 1 ■ P 2 — 2p2'52 x ft ■ pi — 6n*n • P 2 S 2 x n ■ pi — 3re*n • P 1 S 2 x n ■ p 2 


+ 6n*n • P2S2 X n-p2 - 5n*S'2 x pi ■ P2] + 


35G^mi ^ 6G'^m2 ^k 

eikjU *2 4 ^ikj 


n-^Si 


+ 


Gm 2 


[lQeijkp{Sin ■ Pi + beijkU^Sipl - 20eijfcP2>S'i^ • P 2 


Sniir"^ 

+ 2Aeijkn^S\n ■ pin ■ p 2 — x n ■ pi — 24n*n • P 2 S 1 x ft ■ pi + 8 P 2 S 1 x n ■ p 2 


+ 24n*n • P 2 S 1 xn-p 2 + 16n*5i x pi ■ p 2 - 6 eijkn^Si{n ■ ^ 2 )^] + 


7G^m2 
2r3 




H- ^eijkn^Sf, 


mir 


(4.10) 


where we have used the abbreviation pa = Pat'kkia- This results in a compact Hamiltonian: 
Sixn-pi {Apjpi ■ p 2 + IQpi ■ P 2 pI - 16(pi • P 2 f + ^p\ - #2 


uNNLO _Gm2 

"so 
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+ 24pi • np\p 2 ■ n + 2Api ■ np 2 ■ npi ■ p 2 - 12pl{pi ■ n)^ + 3pl{p2 ■ n)^ 

- A8pi ■ P2{P2 ■ n)^ + 2Apl{p2 ■ n)^ + QO{pi ■ n)^iP2 ■ n)^ - 15(p2 • n)"^) 

- 4Si X n • p 2 {2pIpi • P 2 + p\pI + 2pi ■ p 2 pl - 2(pi • ^ 2 )^ + 6pi • np\p 2 ■ n 
+ Qpi ■ np 2 ■ np 2 - 3pl{pi ■ n)^ - 3pl{p2 ■ n)^ + 15(pi • n)‘^{p 2 ■ n)^) 

- 2Si X pi- p 2 {dpi ■ npl - bp\p 2 ■ n + I2pi ■ npi ■ p 2 - lQp 2 ■ npi ■ p 2 


- Api ■ npi + 4p2 • npi + I 2 pi ■ n{p 2 ■ nf - I 8 {p 2 ■ n)^) 

G'^m m r 

- [s. x « ■ A {25p? + 22ft ■ h- 45pl + 76#, • • n - MO® . Hf 

+ 66 (p 2 • n)'^) — 4S'i X n ■ p 2 {llpl — 14pi • p 2 — 16pi • np 2 ■ n — 8(pi • n)^) 


+ 5i X Pi ■ p2 (l65pi • n — 119 p 2 ■ nj 


+ 


G'^m? ■■ 
16r3 . 


Si xn-pi {A2p\ + lOpi • p2 - 57p2 


+ 330pi ■ np 2 ■ n + 48(pi • fi)^ - 229(p2 • n)'^) - Si x n ■ p 2 (l3pf + 228pi • p 2 
— 201^2 + 18pi • np 2 ■ n — 32{pi ■ n)^ + 285(p2 • — 4Si x pi ■ p2{31pi ■ ft 


- 6p2 ■ n) 


+ 


G^niiml 


8 r^ 


1935i X n ■ Pi — 3195i x n ■ p2 


+ 


G^m\m 2 


— 965i X n ■ p2 


• 

9G^mP, 

—♦ —# —* —# ' 

P2 

- 1 - ^ 
47^4 

SS*! X re • Pi — 85i X re • p 2 


8 r^ 

+ (lo2) 


475i X re • Pi 
(4.11) 


The Poisson brackets are the standard canonical ones as shown in [18]. 


4.1 Resolution via canonical transformations 

If the EFT Hamiltonian obtained here in eq. (4.11) is physically equivalent to that of [12, 13], 
then there exists an infinitesimal generator p of a canonical transformation such that 

AH = {H,g} = {Hn + -ffipN + hIq, Pnnlo + 5nlo + 52Pn} 

= AH 2 PN + AHspn + ARNLO ^ ^ ^2) 

where here we have dropped contributions to sectors beyond linear in spin and beyond NNLO, 
and where 

AH = Heft - 77adm- (4-13) 

Thus, the contribution to the NNLO spin-orbit sector comprises 

AH^q^^ = {HvjS'NNLo} + + {-f^LO)5'2PN}5 (4-14) 

and we also have here contributions to lower orders, given by 

AHs^lo = {Hn,4?o}, 

AH 2 PN = {-H'n,52Pn}> (4-15) 

so we also require that the canonical transformation is consistent with the equivalence at 
NLO of the spin-orbit, and 2PN non spinning Hamiltonians. 

Similarly to the construction considerations in [26, 31], we find for the infinitesimal gen¬ 
erator of PN canonical transformations for the NNLO spin-orbit sector, S'nnlO’ following 
general form: 

Pi X P2 (^ 9 ipI + 92P1 ■ P2 + gsPl + 94 (pi ■ + 95 pi ■np2-n 


ffNNLO 


Gm 2 


Si 
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+56 (p 2 -fi^ ^ + Pi X n (^grpjpi ■ n + gsp\p 2 ■ n + g^Pi ■ P 2 P 1 ■ n 
+ gioPi ■ P 2 P 2 ■ n + 511 P 2 P 1 • ^ + 9 i2pIp2 ■ n + gispi ■n[p 2 - 
+514 (pi • n) P 2 -n + gi 5 + pie (p2 ■ ^ + p 2 x n (^gnplpi ■ n 

+ 518 P 1 P 2 • n + P 19 P 1 • P 2 P 1 • n + g 2 oPi ■ P 2 P 2 ■ n + g 2 iplpi ■ n + g 22 plp 2 ■ n 


+523P1 • n I p 2 • n I + 524 fpi • n) P 2 • n + 525 ipi ■ n] + P 26 (p 2 ■ n 


+ 


G^m2 




Pi X P2 (527l«l +528"l2) 


+ Pixn { pi-n { g 29 mi + 5301112) + P 2 ■ n (5311111 + 5321112 
+P 2 X n (Pi • n ( 533 ml + 534m2) + P 2 ■ n (5351111 + 536iii2 


(4.16) 


We should also have the generators contributing first to lower orders, as noted in 
eq. (4.14), so that their coefficients are already set from eq. (4.15). For the NLO spin- 
orbit sector we have the generator from eq. (7.8) in [26, 31] with its coefficients set to the 
values 

51 = 52 = 0, 53 = 54 = 0, 55 = 0. (4.17) 

We should also take into account the generator, which contributes first at the 2PN non 
spinning sector, from eq. (7.10) in [26] with its coefficients set to 

51 =0, 52 = 53 = 0, 54 = 0, 55 = 0, 56 = 0, 57 = (4.18) 

Thus, we plug in eq. (4.14) our ansatz for 5 nnlo (4-16), together with the fixed 

generators in eqs. (4.17), (4.18), and we compare that to eq. (4.13). Comparing 0{G) terms 
fixes the 0{G) coefficients of 5 nnlo values 
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(4.19) 


This eliminates all of the 0{G) terms in the difference. Comparing the remaining O(G^) 
terms in the difference hxes the 0 {G‘^) coefficients of 5 nnlo values 
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(4.20) 

532 - “Y, 

533 = 2, 

534 — “2) 

535 — 0, 536 — “^• 

This eliminates all of the O(G^) 

terms, as well as 

all terms at 0{G^) in the difference. 

Hence, 


we have shown that the ADM Hamiltonian and the EFT Potential at NNLO spin-orbit are 
completely equivalent. 
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4.2 Hamiltonians in the center of mass frame 


Since the center of mass of a binary moves uniformly along a straight line, it makes sense 
to separate this motion and reduce the number of orbital variables. This is achieved by a 
specialization of the Hamiltonian to the center of mass frame, where the total linear momen¬ 
tum vanishes. Here the total linear momentum is just the sum of the individual momenta, 
P 1 +P 2 = 0, since we are considering the conservative sector where the recoil due to gravita¬ 
tional waves vanishes. The orbital dynamics is then governed by a single canonical momentum 
p = Pi = —p2 and its conjugate r = yi — y2, the binary separation vector. 

We further introduce the tetrad basis h, A, I, where I = L/L, X = I x ft, and the orbital 
angular momentum is L = rn x p, where a similar tetrad was used in [38] for an approximate 
solution of the EOM with spin. This is a sensible choice for the description of a spinning 
binary, as it allows, for instance, a separation of p into circular and radial components. 


p = PrU -I—A, 
r 

where pr = p ■ n. However, we will not expand the spin products as in 
Sa- Sb = SA-nn- Sb + Sa-XX-Sb + Sa- U- Sb, 


(4.21) 


(4.22) 


since these spin products better highlight the structure of the interaction and keep expres¬ 
sions more compact. Furthermore, the products Sa ■ Sa are constants in this work. As a 
convention, we then replace multiple products involving L by the products Sa - Sb, so with 
this prescription at most one scalar product containing L remains in each term. 

As in [26] we are expressing the results in terms of dimensionless variables, where di¬ 
mensions of length are rescaled by the total mass Gm = G{mi +m 2 ) and dimensions of mass 
are rescaled by the reduced mass p = mim 2 /m. This rescaling is denoted by a tilde, e.g., 
L = The dependence on the masses is then reduced to the mass ratio q = mi/m 2 , or 

the symmetric mass ratio v = p/m. 

We use here the Hamiltonians arising from the formalism recently set up in [18], and 
hence within the corresponding gauge choices. We include also the results from the cubic 
order in spin from [32], see also [39, 40] for the black hole case, and [41]. All expressions 
here are then complete to the 3.5PN order with spins. Before we present the spin-dependent 
Hamiltonians, we provide for completeness the point-mass Hamiltonians to 2PN order in the 
same formalism [35] and center of mass frame, reading 
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(4.25) 


The 2PN order in the point-mass sector is sufficient for deriving the spin-dependent part of 
the binding energy to 4PN order, or to transform the spin-dependent Hamiltonians to FOB 
gauge up to 4PN order. 
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The center of mass spin-orbit Hamiltonians to 3.5PN order read 
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where the NNLO Hamiltonian was derived in the present paper and the lower order Hamil¬ 
tonians are taken from [18]. The relevant Hamiltonians at quadratic order in spin are also 
given in [18]. After transformation to the center of mass form, they read 
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(4.32) 


The LO cubic in spin Hamiltonian at 3.5PN order from [32] reads in the center of mass frame 
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Recall that in the black hole case it holds = 1 = ^^0 that [1 o 2] implies 

q o q~^. It is nice to note that in the representation that we picked, the Hamiltonian H^q 
appears as a correction to the spin-orbit interaction at 3.5PN order, which is nonlinear in the 
spins. It should also be noted though that the representation is not unique. One can expand 
a scalar product of spins with eq. (4.22), which produces two products involving L, and then 
absorb two different products involving L by rewriting them in terms of a spin product. 

All of the results in the present section are still gauge dependent. We are going to work 
out gauge invariant relations in the following section, which however require restrictions on 
the orbit and spin orientations. 


5 Complete gauge invariant relations to 3.5PN order with spins 


Starting from the Hamiltonians given in the last section, it is straightforward to obtain gauge 
invariant relations between the binding energy and the orbital frequency or orbital angular 
momentum for circular orbits. Circular orbits are defined by f = const. This holds if the 
radial momentum vanishes, p,- = 0. Since there is no obvious way to define spin orientations 
in a gauge invariant manner, we further restrict to aligned spins, Sa ■ n = ^ = Sa ■ X, and 
Sa ■ L = SaL. The relations presented here were given in part already in [26]. Here we 
complete these relations to the 3.5PN order with spins by including recent results for the 
cubic order in spin [32]. The computation follows that in [26]. 

The circular orbit relation between L and f takes the form 
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(5.1) 


This relation is gauge dependent and therefore differs from the result in ADM gauge given by 
eq. (8.22) in [26]. A gauge invariant quantity is the orbital frequency lo, which is commonly 
expressed in terms of x = in the binding energy. The gauge invariant relation between 
L and x has the following addition compared to eq. (8.24) in [26]: 
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Finally, the gauge invariant relation for the binding energy as a function of the orbital 
angular momentum reads 
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and as a function of x it reads 
espin(3:) = vx^!"^ (^1 + 52 
+ vx 


4 

- 3 +" 


31z^ , \ 2 ( 

—-41-.U 


/7i/2 211i2 27 


(^5i/g + § 2 ( 1 ^ 


-i+Af-f 


V 12 

2 


+ T 


3912 27 


+ 12 X^SiS2 


1 + X 


5z2 5 

18 6 


+ <[ 12 x ^82 


V 8 

7s5!) + ^ _ 3) 


6 


+ 


•'2(ES2) 


{12 + 1) ) ( Y^ + 


25 5C, 


2(ES2) 




^52 [2(<72(es2) - <^2(BS3)) 


+ g(3<72(E32 ) 2(72(3g3))] + [(?(6 4(72(Eg2)) + (4 3(72(Eg2))] 


+ [1^2] 


(5.4) 


The point-mass gauge invariant relations to 4PN order can be found in [6]. 

6 Conclusions 

In this work we implemented the EFT for gravitating spinning objects in the PN scheme 
[18] at the NNLO level, which was first treated in [19]. We derived the NNLO spin-orbit 
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interaction potential at the 3.5PN order for rapidly rotating compact objects. Such high 
PN orders are required for the successful detection of gravitational radiation, as the EOB 
Hamiltonian, e.g., requires parameters for the even higher 5PN and 6PN orders in the point- 
mass case, in order to produce good waveforms. From the NNLO spin-orbit interaction 
potential, which we obtain here in a Lagrangian form for the first time, we directly derived 
the corresponding Hamiltonian. We then compared our result to the ADM Hamiltonian result 
[12], and arrived at a complete agreement between the ADM and EFT results. Therefore, in 
order to complete the spin dependent conservative sector to 4PN order, it remains to apply 
the EFT for gravitating spinning objects [18] at NNLO to quadratic level in the spin, and for 
higher order in spin finite size effects, as was already done in [32] for cubic and quartic orders 
in spin. The NNLO spin-squared result is indeed presented in another recent paper [27], 
which then completes the conservative sector to 4PN order. Finally, we provide the relevant 
Hamiltonians in the center of mass frame, and the complete gauge invariant relations among 
the binding energy, angular momentum, and orbital frequency of an inspiralling binary with 
generic compact spinning components to 3.5PN order. 

The spin-orbit sector constitutes the most elaborate spin dependent sector at each or¬ 
der, and accordingly we encountered here a proliferation of the relevant Feynman diagrams, 
where there are 132 diagrams contributing to this sector, and a significant increase of the 
computational complexity, e.g. there are 32 two-loop diagrams here. We also recall that as 
the spin is derivative-coupled, higher-order tensor expressions are required for all integrals 
involved in the calculations, compared to the non spinning case. However, the computa¬ 
tion is made efficient through the use of the “NRG” fields, which are advantageous also in 
the spin dependent sectors, together with the various gauge choices included in the EFT 
for gravitating spinning objects [18]. In addition, we automatized the EFT computations 
here, and carried out the automated computations in parallel. Hence, it is clear that for 
higher order corrections automated EFT computations, utilizing the “NRG” fields, should 
be implemented, and are most powerful and efficient. It should be stressed that in order to 
obtain such higher order results, all lower order results are required consistently within one 
formalism, and so also for that the derivation presented in this work is essential. This work 
then paves the way for the obtainment of the next-to-NNLO spin-orbit interaction potential 
at 4.5PN order for rapidly rotating compact objects, once this level of accuracy would be 
approached. 
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A Irreducible two-loop tensor integrals 

In the evaluation of the irreducible two-loop diagrams we encounter irreducible two-loop 
tensor integrals up to order 3. These are reduced using the integration by parts method to a 
sum of factorizable and nested two-loop integrals, as explained in [19], and see appendix A 
there. In addition to the irreducible two-loop tensor integrals, which were given in appendix A 
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of [19], the two following reductions are also required here: 
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